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■ ABSTRACT 



The first observations conducted as part of the Chandra ACIS survey of M 33 (ChASeM33) 
sampled the eclipsing X-ray binary M 33 X— 7 over a large part of the 3.45 d orbital pe- 
riod and have resolved eclipse ingress and egress for the first time. The occurrence of the X- 
ray eclipse allows us to determine an improved ephemeris of mid-eclipse and binary period as 
HJD (245 3639.119 ± 0.005) ±N x (3.453014 ± 0.000020) and constrain the eclipse half angle to 
26.5° ±1.1°. There are indications for a shortening of the orbital period. The X-ray spectrum is 
best described by a disk blackbody spectrum typical for black hole X-ray binaries in the Galaxy. 
We find a flat power density spectrum and no significant regular pulsations were found in the 
frequency range of 10~ 4 -0.15 Hz. HST WFPC2 images resolve the optical counterpart, which 
can be identified as an 06III star with the help of extinction and colour corrections derived from 
the X-ray absorption. Based on the optical light curve, the mass of the compact object in the 
system most likely exceeds 9Mq. This mass, the shape of the X-ray spectrum and the short term 
X-ray time variability identify M 33 X— 7 as the first eclipsing black hole high mass X-ray binary. 

Subject headings: Galaxies: individual: M 33 — X-rays: individuals: M 33 X— 7 — X-rays: binaries — 
binaries: eclipsing 
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1. Introduction 

M 33 X— 7 (hereafter X— 7) was detected as a 
variable source by the Einstein observatory with 
a maximum luminosity in the 0.15-4.5 keV band 
(assuming an absorption column of 10 21 cm~ 2 ) 
that exceeds 10 38 erg s _1 (Long et al. 1981; Mark- 
ert & Rallis 1983). The source stayed active 
in all following observations. Its variability was 
explained by an eclipsing X-ray binary (XRB) 
with an orbital period of 1.7 d and an eclipse 
duration of -0.4 d (Peres et al. 1989; Schul- 
man et al. 1993, 1994). Based on ROSAT and 
ASCA data (Larson & Schulman 1997; Dubus 
et al. 1999, hereafter DCL99), the orbital pe- 
riod was found to be twice as long. DCL99 de- 
scribed the shape of the eclipse by a slow ingress 
(A$i, lgross = 0.10 ± 0.05), an eclipse duration of 
A$ccii P sc = 0.20 ± 0.03 and a fast eclipse egress 
(A$ ogross = 0.01 ± 0.01) with an ephemeris for 
the mid-eclipse time of HJD 2448631. 5±0.1 + 
Nx(3.4535±0.0005). DCL99 also report 3cr ev- 
idence for a 0.31 s pulse period. They come 
to that conclusion by splitting PSPC and HRI 
data of 500 s intervals of continuous data where 
X— 7 was positively detected and calculating the 
summed Rayleigh power spectrum. The linearly 
binned power spectrum showed an significant ex- 
cess (at 99.9 per cent confidence) at the proposed 
period when compared to simulated data assum- 
ing a Poisson distribution with the X— 7 mean flux. 
The signal is broader than what would be expected 
from a simple sinusoidal pulse. It was not possi- 
ble to check if the variability is coherent. DCL99 
conclude "Although pulsed emission seems a rea- 
sonable assumption, the power excess could arise 
from variability of a different nature (e.g. broad- 
band variability increasing the chance of spurious 
detection)." The orbital period, pulse period and 
observed X-ray luminosity are remarkably similar 
to those of the Small Magellanic Cloud neutron 
star XRB SMC X-l (Liu et al. 2000). X-7 was 
the first and only identified eclipsing accreting bi- 
nary system with an X-ray source in an external 
galaxy other than the Magellanic Clouds before 
the detection of similar behavior based on XMM- 
Newton and Chandra data of the NGC 253 X-ray 
source RX J004717.4-251811 (Pictsch et al. 2003). 

Pietsch et al. (2004, hereafter PMM2004) an- 
alyzed several observations of the XMM-Newton 



M 33 survey and an archival Chandra observation 
where X— 7 was in the field of view. The observa- 
tions cover a large part of the 3.45 d orbital period, 
however not eclipse in- and egress. PMM2004 de- 
tected emission of X— 7 during eclipse and a soft X- 
ray spectrum of the source out-of eclipse that can 
best be described by brcmsstrahlung or disk black- 
body models. No significant regular pulsations of 
the source in the range 0.25-1000 s were found. 
The average source luminosity out of eclipse was 
5xl0 37 erg s^ 1 (0.5-4.5 keV, corrected for Galac- 
tic foreground absorption). In a special analysis 
of DIRECT 1 observations PMM2004 identified as 
the optical counterpart a B0I to 071 star of 18.89 
mag in V which shows the ellipsoidal light curve 
of a high mass X-ray binary (HMXB) with the 
X— 7 binary period. Based on the location of the 
X-ray eclipse and the optical minima, PMM2004 
derived an improved binary ephemeris and argued 
that the compact object in the system is a black 
hole. Those authors reached this conclusion based 
on the mass of the compact object derived from or- 
bital parameters and the optical companion mass, 
the lack of pulsations, and the X-ray spectrum. 
X— 7 would be the first detected eclipsing high 
mass black hole XRB. 

The Chandra ACIS-I survey of M 33 (ChASeM33) 
is a very large program which will accumulate in 
seven deep pointings, each 200 ks in length, a to- 
tal exposure of 1.4 Ms. During several of these 
pointings X— 7 is in the field of view. We report 
here on ChASeM33 observations of X— 7 spread 
over just 20 binary orbits, which resolved for the 
first time the eclipse ingress and egress and al- 
lowed us to constrain the light curve of X— 7 for 
binary phases around eclipse. Preliminary results 
of the first observations were announced by Sasaki 
et al. (2005). In addition we identify the source 
on archival HST WFPC2 images. Throughout the 
paper, we assume a distance to M 33 of 795 kpc 
(van den Bergh 1991). 

2. Chandra observations and results 

X— 7 was sampled by the ACIS-S3 chip in one 
Chandra observation and by the 10, II and 12 
chips during 11 additional observations. Table 1 
summarizes these observations giving observation 

1 For information on the DIRECT project see 
http : //cf a-www . harvard . edu/~kstanek/DIRECT/ . 
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identification (ObsID) in column 1, observation 
start date (2), elapsed time (3), the ACIS chip cov- 
ering X— 7 (4), the offset of X— 7 from the pointing 
direction (5), and the X— 7 binary phase during 
the observation (6) using ephemeris that will be 
discussed in Sect. 2.1 and 4.1. The source bright- 
ness varied from ^3xl0~ 3 ct s _1 to ^0.2 ct s -1 
normalized to ACIS-I on-axis. 

Of the pointings considered here, the Field 5 
position of X— 7 is located close enough to the opti- 
cal axis to result in significant pile-up, particularly 
during the high phase. To correct the light curve 
for the effects of pile- up, we used the best-fit model 
for unpiled data from the other pointings (see 
Sect. 2.2), and applied a pileup(phabs(diskbb)) 
model within XSPEC. We estimated the pile- 
up parameters by fitting simultaneously the high 
state data for ObsID 6382 (insignificant pile-up) 
and ObsIDs 7170 and 7171 (significant pile-up). 
The phabs and diskbb models were set to their 
best fit values from the unpiled spectra; the diskbb 
normalizations were tied for all three datasets. All 
of the pile-up model parameters except the grade 
morphing parameter a were frozen. For a descrip- 
tion of the pile-up model see Davis (2001). The 
pile-up parameters were set to their default val- 
ues, but the fr_time parameter was frozen at 3.2 
for ObsIDs 7070, 7071, and at for 6382; this has 
the effect of turning off the pile-up model compo- 
nent for ObsID 6382. The resulting reduced \ 2 of 
the fit was 0.97 and the corresponding value for 
a was 0.692t°;^ (90% confidence limits). Using 
the fitted value for a, a correction curve was gener- 
ated by evaluating the ratio of the XSPEC "Model 
predicted rate" without pile-up (fr_time=0) to the 
corresponding rate with pile- up (fr_time=3.2) as 
a function of the model predicted rate (with pile- 
up). This ratio was used to correct the observed 
count rates of ObsIDs 6384, 7070, and 7071 for the 
effects of pile-up. 

For the spectral analysis we used standard Level 
2 event files cleaned for columns with higher back- 
ground rates than adjacent columns. In the ob- 
servations of Fields 4 and 5, counts from X— 7 are 
rejected when the source is moving across rejected 
columns (due to the satellite dithering). This ef- 
fect creates spurious periods in period analysis and 
can reduce the counts in 1000 s integration inter- 
vals by varying amounts (up to more than 30%). 
Therefore, we created new Level 2 event files for 



broad band time variability analysis that did not 
reject these columns. We correct for the detection 
efficiency at different off-axis angles using factors 
derived from response files for the different CCDs 
and off-axis positions assuming the disk-blackbody 
spectrum derived in Sect. 2.2. We normalized the 
count rates to a CCD 3 on-axis rate. 

The data analysis was performed using tools in 
the ESO-MIDAS v05SEPpll.O, EXSAS v03OCT_EXP, 
CIAO v3.2 and LHEASOFT v5.3 software pack- 
ages as well as the imaging application DS9 
v3.0b6. 

2.1. Time variability 

During ObsID 6378 and 7171 we observed tran- 
sitions by X— 7 into and out of eclipse, respec- 
tively. We sampled both background and solar 
system barycenter corrected light curves of X— 7 
with a time resolution of 1000 s. To increase the 
signal to noise specifically in the far off-axis Ob- 
sID 6378, we restricted the analysis to the 0.5- 
5 keV band which covers most of the source flux. 
Further sub-dividing this energy band into a hard 
and soft band did not show any significant hard- 
ness ratio changes of the in- and egress behavior 
as could be expected based on EXOSAT observa- 
tions of LMC X-4 and Her X-l (see Figs. 78 and 
79 in Dennerl 1991). To determine eclipse start 
and end times, we approximated the light curves 
assuming constant count rates within and out-of 
eclipse and a linear transition in between using a 
X minimization technique. We searched for la 
errors of the eclipse start and end time, respec- 
tively, assuming those to be the only interesting 
parameter for the fit. 

In ObsID 6378 the transition into eclipse lasted 
12.75 ks and ended at HJD 245 3635.41 10±0.0037. 
The transition out of eclipse in ObsID 7171 
started at HJD 245 3642.8272±0.0052 and lasted 
for 10.52 ks. From these ingress and egress times 
separated by two orbital periods, we directly de- 
rive the mid eclipse ephemeris of the eclipse in 
between to HJD 245 3639.119 ± 0.005. Assum- 
ing this epoch as phase zero and a binary pe- 
riod of 3.453014 d (see Sect. 4.1), we calculated 
light curves of all X— 7 observations (Fig. 1). We 
also determine an eclipse duration of less than 
0.147±0.006 in phase corresponding to an eclipse 
half angle of 26.5° ± 1.1°. 
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ObsIDs 7196 and 7199 fully fall into eclipse. 
ObsID 6384 at phase 0.83 to 0.90 indicates a much 
longer transition into eclipse than ObsID 6378. 
ObsID 7198 shows dipping behavior well before 
eclipse with a return to the out-of eclipse inten- 
sity level at the end. ObsID 6382 shows a sec- 
ond egress from eclipse which is significantly faster 
than the one 16 orbits earlier (ObsID 7171). How- 
ever, the observation starts during egress and the 
phase range within eclipse is not covered. Ob- 
sID 7226 covers eclipse ingress at the end of the 
same orbit and again shows strong dipping well 
before the ingress. Eclipse egress and ingress 
times are consistent with the times derived above. 
Generally speaking the variability of X— 7 before 
eclipse (phase 0.7 to 0.9) seems to be much more 
pronounced in individual observations than after 
eclipse (phase 1.1 to 1.4). Average count rates 
in eclipse and out of eclipse are 0.003 ct s _1 and 
0.15 ct s _1 , respectively. This out-of-eclipse count 
rate varies in different binary orbits by factors of 
1.3 and there are residual short term fluctuations 
that can be described by dips with a similar ampli- 
tude and a duration of several 1000 s also outside 
the pre-eclipse phase. 

To search for pulsations we extracted light 
curves in the 0.5-5 keV band from the longest 
observations after eclipse, ObsIDs 6382 and 7170. 
We created power density spectra in the frequency 
range of 10~ 4 -0.15 Hz and found no significant 
periodic signal with a 3cr upper limit of 5.3% for 
sinusoidal variations. A power density spectrum 
derived from ObsID 6382 by adding the power 
spectra from 23 intervals of 3319 s length (1024 
time bins with the instrument resolution of 3.241 
s) is shown in Fig. 2. The power spectrum is flat 
at a value of 2. 

2.2. Energy spectra 

We analyzed energy spectra of X— 7 for all ob- 
servations. For ObsIDs 6378 and 7171, times 
out of eclipse, during ingress and egress, and 
in eclipse were handled separately. Absorbed 
power-law, brcmsstrahlung and disk-blackbody 
models (which were found to best represent the 
XMM- Newton and Chandra spectra analyzed by 
PMM2004) were first fit to the individual spec- 
tra. The power-law fit to the spectrum obtained 
from ObsID 6382 (which has the highest statis- 
tical quality) yields an unacceptable fit with a 



reduced % 2 of 1.92 while bremsstrahlung and disk- 
blackbody models result in \r °f 1-45 and 1.32, 
respectively. The derived parameters are consis- 
tent within the errors for all spectra including the 
spectra accumulated during eclipse ingress and 
egress and during ObsID 7198 when the source 
shows high variability. Therefore, we performed 
a simultaneous fit with the disk-blackbody model 
to the spectra of eight out-of-cclipse observations 
excluding times when the X-ray source was in 
eclipse and forcing the absorbing column density 
to be the same for all observations (this reduces 
the number of free fit parameters compared to the 
individual fits). The resulting inner disk temper- 
ature was systematically higher (at ^1.3 keV) for 
ObsIDs 6384, 7170 and 7171: these are the obser- 
vations which include X— 7 nearly on-axis. The 
most likely reason for the "harder" spectra is pile- 
up and we therefore excluded these three spectra 
from further analysis. The fit to the remaining five 
spectra from ObsIDs 6378, 6382, 6386, 7197 and 
7198 shows no significant differences in the inner 
disk temperature and therefore this parameter was 
also forced to be the same in the simultaneous fit. 
We refit the power-law and bremsstrahlung mod- 
els for comparison. The best fit is obtained with 
the disk-blackbody model yielding a \r °f 1-10. 
The Xr f° r the best fitting bremsstrahlung and 
power-law models are 1.16 and 1.44, respectively. 
The derived spectral parameters are similar to the 
ones reported by PMM2004: disk-blackbody with 
inner disk temperature kT = 0.99±0.03 keV and 
N ff = (0.95±0.10)xl0 21 cm" 2 ; bremsstrahlung 
with temperature kT = 2.74±0.13 keV and N ff = 
(2.05±0.12)xl0 21 cm~ 2 ; power-law with photon 
index 7 = 2.38±0.05 and N ff = (3.32±0.17) x 10 21 
cm~ 2 . The best fit disk-blackbody model is shown 
in Fig. 3. 

The normalization of the disk-blackbody model 
is given by K = (r-m/d) (cos i) with the inner disk 
radius r; n in km, the source distance d in units 
of 10 kpc and the disk inclination i. The spectra 
show variations by a factor of ~2 in normalization 
with K=0.054 for ObsID 6378 and relative factors 
0.58, 0.94, 1.13 and 1.37 for ObsIDs 7198, 7197, 
6386 and 6382, respectively. 

The N# values for the models discussed above 
clearly indicate absorption within M 33 or intrin- 
sic to the source in addition to the Galactic value 
(5.86 and 6.37 xlO 20 cm" 2 in the direction ofX-7 
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according to Dickey & Lockman 1990; Stark et 
al. 1992, respectively). Absorbed and unabsorbed 
source fluxes in the 0.3-10 keV band are in the 
range (5.4-12.6) x 10 -13 erg cm~ 2 s _1 and (6.2- 
14.7) xl0~ 13 erg cm~ 2 s _1 , respectively, based on 
the best fitting disk-blackbody model. These 
fluxes correspond to source luminosities of (4.1- 
9.6)xl0 37 ergs" 1 and (4.7-11.2) x 10 37 erg s" 1 , 
respectively. 

The N/f value of the best fitting disk-blackbody 
model indicates that X— 7 lies on the near side of 
M 33 as the absorbing column within M 33 can 
be determined to ^2.2 xlO 21 cm -2 from a 47 x 93 
arcsec half power beam width Hi map (Newton 
1980). From the value we can compute the 
expected optical extinction Ay — 0.53 ± 0.06 mag 
and E(B-V) = 0.18 ±0.02 using the standard re- 
lations (Predehl & Schmitt 1995). These numbers 
are in the range given by PMM2004 who assumed 
that we see X— 7 through less than half the ab- 
sorbing column within M 33. 

2.3. Improved position 

X— 7 was located in four of the five ChASeM33 
fields observed so far. The source was closest to 
on-axis (1.58' off-axis) in the Field 5 observations 
(see Table 1). This results in the most compact 
PSF and thus the most reliable position deter- 
mination. To refine the absolute astrometry of 
the Field 5 data, we searched the USNO-B1.0 and 
2MASS catalogs for close positional matches with 
X-ray sources. Similarly, to improve the statistics 
for X-ray ccntroiding, we worked with a merged 
dataset for ObsIDs 6384, 7170, and 7171. We 
identified 9 candidate optical/2MASS objects. Six 
were rejected because of far off-axis positions or 
small number of counts. The remaining candi- 
dates were > 7' off-axis, for a 2MASS object (~250 
counts, 0.5-5 keV) which was 7.8' off-axis. We 
enhanced the number of candidates for registra- 
tion by adding two isolated centrally brightened 
supernova remnants (SNRs) with good radio posi- 
tions (sources 57 and 64 from the list of Gordon et 
al. 1999) assuming that the finite size of SNRs - 
and potential differences in the X-ray versus radio 
distribution - do not bias the position determina- 
tion. We determined the X-ray centroids based 
on an iterative sigma-clipping algorithm applied 
to the 0.5-5 keV X-ray data. Based on an initial 
position estimate and clipping radius, the stan- 



dard deviation of the radial distribution is evalu- 
ated, and points greater than a given number of 
standard deviations are rejected. The iteration 
of centroiding and rejecting events continues until 
the centroid converges to within a specified tol- 
erance or for a fixed number of iterations (10). 
The difference in sky coordinates between the cat- 
alog position and the X-ray centroid position was 
evaluated for each source. The mean offset was 
Ax = (0.23±0.37)", Ay = (0.32±0.21)". The cen- 
troid of the X— 7 source was evaluated in the same 
way, and the resulting offset was applied to correct 
the sky position. Finally, the corresponding ce- 
lestial coordinates were evaluated using the CIAO 
tool dmcoords; the aspect solution files (asoll files) 
were used in order to correct the positions using 
the aspect offsets. The resulting X— 7 position 
is RA.J2000 = 01 /l 33 m 34!12, (Jjaooo = +30°32'11 , .'6 
with a combined error of 0.5". The position is 
within 2cr of that given by PMM2004 based on a 
registration using just one SNR. 

3. Optical observations and results 

M 33 X— 7 is located along the line of sight 
to the dense OB association HS 13 (Humphreys 
& Sandage 1980): it was identified by PMM2004 
with a specific star within this association, as a re- 
sult of detection of regular (ellipsoidal) variations 
in B and V at the X-ray period. The OB asso- 
ciation has been imaged with HST using WFPC2 
in three filters: F336W, F439W, and F555W. The 
observing details are listed in Table 2. In an at- 
tempt to learn more about the optical counterpart 
to X— 7, we retrieved relevant datasets from the 
MAST archive. Data retrieved from the archive 
are automatically reprocessed with the latest cali- 
bration files. After inspecting the individual expo- 
sures, we combined the two exposures taken with 
the F336W filter and the two taken in F439W (us- 
ing the STSDAS task gcombine) to eliminate the 
effects of cosmic rays on the images. 

The position of the X-ray source as determined 
with Chandra is shown in Fig. 5. Since the er- 
ror in the astrometric solution of images in the 
HST pipeline is typically <~1.5"— 2", we attempted 
to reduce the positional uncertainty by register- 
ing the WFPC2 images to the USNO-B1.0 or 
2MASS frames. Given the small field of view and 
greatly superior spatial resolution of the WFPC2 
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images, many sources identified in the USNO-B1.0 
or 2MASS catalogues are resolved into multiple 
objects, making a unique match between stars in 
the image and sources in the catalogue difficult. 
Therefore, we performed our astrometric correc- 
tions in two steps, taking advantage of the much 
larger field of view of the KPNO Mosaic B image 
of M 33 (Massey et al. 2002). First, we identified 
a sample of 25 isolated, bright USNO-B1.0 stars 
within 4.'5 of X— 7 to compute the shift required 
to bring the Mosaic B image into the USNO-B1.0 
frame. The RMS positional error of these stars 
in the corrected Mosaic image was unacceptably 
large (0'.'85). We then tried the same procedure 
using 21 bright 2MASS stars located in isolated 
areas in the KPNO Mosaic B image. The RMS 
positional error of these stars in the corrected Mo- 
saic image was 0'.'16. We then used 10 bright, 
isolated stars in common between the Mosaic B 
and F439W images (restraining the selection to 
the region covered by the WFPC2 CCD in which 
the X— 7 counterpart was located) to compute the 
shift required. The RMS positional error in the 
shifted HST image was 0'.'13. Combining the 0"1 
absolute uncertainty of the 2MASS position with 
the uncertainties listed above, the absolute astro- 
metric uncertainty of the final registered F439W 
image is 0'.'23. 

We applied the calculated shifts above to the 
HST image. Fig. 5 shows a 10" x 10" field from 
the F439W image centered on M 33 X-7. The 
black solid circle shows our best estimated posi- 
tion for the X-ray source from the Chandra image. 
The error circle of X— 7 is 0'.'5 in radius (see above). 
The HST positional accuracy is indicated by a 
cross. We also show the error quoted by PMM2004 
in red. The error circles overlap with a bright 
star, coincident with that proposed by PMM2004 
(based also on time variation arguments) as the 
donor star for the compact object. 

We also carried out aperture photometry of the 
stars in the HST field (using the IRAF procedures 
daofind and phot). We find that the optical coun- 
terpart to X— 7 with the WFPC2 resolution is not 
a blend of stars. PSF fits to the source in the 
F336W, F439W, and F555W images give FWHM 
compatible with the other point-like sources in the 
images. The star about 0"9 to the South on the 
other hand is just resolved in at least two equally 
bright sources which are separated by ~0f.'2. We 



can rule out the presence of another star of similar 
brightness at the position of the optical counter- 
part of X— 7 to this distance, which corresponds 
to a projected separation at the distance of M 33 
of 0.8 pc. 

The optical counterpart has apparent magni- 
tudes of 17.6, 18.2, and 18.9 for the F336W, 
F439W, and F555W filters, respectively, in the 
STMAG system. The colors derived from these 
filters for the counterpart to X— 7 are typical of 
the other bright stars (my < 20.5) in HS 13 as 
observed with WFPC2. The F336W, F439W, and 
F555W filters are centered approximately on the 
corresponding U, B and V filters. Adopting a color 
transformation of 0.5 mag for U, 0.66 mag for B, 
and 0.03 mag for V (Mobasher 2002), we find mu 
of 18.1 mag, me of 18.8 mag, and my of 18.9 mag 
(i.e. U-B of -0.7 mag, B-V of -0.1 mag) for the 
optical counterpart to X— 7. 

4. Discussion 

The well sampled orbital light curve of X— 7 
indicates stronger variability before eclipse com- 
pared to after eclipse (Fig. 1). Variability at this 
phase is often observed in HMXBs and is explained 
by the viewing geometry through the innermost 
regions of the wind of the companion and dense 
material following the compact object in its or- 
bit (e.g. Haberl & Day 1992). Dense structures 
are created by the gravitational and radiative in- 
teractions of the compact object with the stellar 
wind (Blondin et al. 1990, 1991). This behavior 
is also reflected in the on-average longer eclipse 
ingress time and longer eclipse duration derived 
by DCL99. 

Similarly to PMM2004 we find residual emis- 
sion from the source during eclipse. Residual emis- 
sion during eclipse was measured from most eclips- 
ing XRBs and can be explained by re-processing of 
primary photons from the compact X-ray source 
in an extended accretion disk corona (which is not 
fully occulted) or by scattering in the compan- 
ion atmosphere/stellar wind. Residual emission 
of up to ~10% of the uneclipsed flux was reported 
(Haberl 1991; Lewis et al. 1992; Ebisawa et al. 
1996) depending on system geometry and wind 
density. The X— 7 residual emission of ~4% is well 
within these limits. 
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4.1. Improved ephemeris 

DCL99 modeled the folded light curve from 
X— 7 as a constant flux plus linear ingress and 
egress plus an eclipse interval with zero flux. It 
is obvious from their data (see their Fig. 1) that 
the eclipse egress is better determined than the 
eclipse center and duration. These strongly de- 
pend on the shape of the pre-eclipse dips contained 
in the light curve as can be seen from the resolved 
pre-eclipse behavior in the Chandra data (Fig. 1). 
DCL99 did not determine the eclipse parameters 
from individual eclipses but only from the aver- 
age light curve due to limited statistics. The time 
of eclipse egress is the best determined parame- 
ter. Unfortunately, in the paper they do not give 
this parameter separately but only the center of 
eclipse and length of eclipse. In the following we 
use eclipse egress times to determine an improved 
orbital period P and a possible period derivative 
P. 

Due to limited phase coverage, the XMM- 
Newton observations do not resolve eclipse ingress 
or egress. The time of eclipse egress can only 
be constrained to <0.02 in phase. Based on 
these data PMM2004 restricted the time of eclipse 
egress to HJD 245 1760.953 ± 0.035 (note typo- 
graphical error in egress time in PMM2004, but 
the calculation used the number given here) and 
determined a time of mid-eclipse assuming the 
eclipse shape parameters of DCL99. With this 
mid-eclipse epoch and the one given by DCL99, 
PMM2004 determined an improved orbital period. 

The maximum eclipse duration of 0.147±0.006 
(Sect. 2.1) determined from individual observa- 
tions is - as expected - significantly shorter than 
the one given by DCL99 (0.20±0.03) from average 
orbit fitting. From the DCL99 ephemeris it is not 
possible to derive a well defined eclipse egress time. 
The same is true for the Einstein results (Peres et 
al. 1989, hereafter PRC89). We therefore decided 
to re-analyze relevant ROSAT and Einstein data. 

ROSAT PSPC observations of M 33 were short 
compared to the HRI observations and did not 
cover the X— 7 eclipse egress. We therefore re- 
stricted our re-analysis to ROSAT HRI data. 
After screening for high background, we com- 
bined data with continuous observation intervals 
which lead to variable integration times of typ- 
ically 1800 s (minimum 41 s, maximum 3791 s) 



depending on the duration of the scheduled ob- 
servation and background. DCL99 in contrast 
grouped the data in 3000 s averages. Only once 
is the eclipse egress closely monitored. In the 
ROSAT interval corresponding to ObsID 600488h, 
X-7 was still in eclipse while in the following inter- 
val (corresponding to ObsID 600489h), less than 
0.04d later, X-7 already featured a count rate that 
indicated the source was out of eclipse. This al- 
lowed us to restrict the time of eclipse egress to 
HJD 244 9571. 724± 0.018. In another case during 
ObsID 600020h-l, observations in eclipse and out 
of eclipse are separated by 0.26 d, which does not 
allow us to further constrain eclipse egress times. 

PRC89 reported Einstein IPC and HRI obser- 
vations of X— 7. For the HRI observations PRC89 
combined several continuous observation intervals 
to get significant data. For the IPC observations 
PRC89 simply integrated over individual contin- 
uous observations. Due to the better statistics 
in the IPC observations we considered only IPC 
eclipse egress coverages (i.e. ObsID 12090, see 
also PRC89, Fig. 1). In images of the observation 
X— 7 is only visible in the energy band 0.6-2.8 keV 
(PI 4-9). We therefore restricted the analysis to 
this energy band. We selected extraction position 
and area by comparison with the close-by bright 
central source X— 8. We specifically investigated 
the "first rising episode at day 1.5" as identified 
by PRC89. In contrast to the report by PRC89, 
X— 7 count rates during this period do not show 
increasing flux but the source intensity is compat- 
ible with zero during all three intervals. In the 
next set of observation intervals at around day 2 
the source is clearly out of eclipse. This indicates 
an eclipse egress between HJD 2444087.840 and 
HJD 2444088.255. 

Combining the ROSAT eclipse egress bound- 
aries with the Chandra ephemeris suggests an or- 
bital period of 3.453014±0.000020 d. The lower 
panel of Fig. 1 shows the ROSAT HRI light curve 
folded with the above period and assuming the 
Chandra mid-eclipse epoch, for the same phase 
range as the Chandra data above. As can be seen 
in Fig. 4, these ephemerides also are consistent 
with the boundaries determined for the XMM- 
Newton eclipse egress. However, they seem to 
miss the boundaries determined above for the Ein- 
stein IPC observation. If we assume a constant 
rate of change of the orbital period over this time, 
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we can model all eclipse egresses. The parabola 
shown in Fig. 4 assumes a period of 3.45294 d 
at the Chandra epoch and orbital period decay 
rate of P or b/-Porb = — 4 x 10~ 6 yr -1 . It nicely 
models the average eclipse egress times of XMM- 
Newton, ROSAT and Einstein. However, also a 
small Porb/^orb = —0.7 x 1CT 6 yr _1 and a pe- 
riod of 3.45302 d or a much higher P rb/-Porb = 
-7.5 x 10~ 6 yr" 1 and a period of 3.45285 d would 
still be consistent with all the eclipse egress bound- 
aries. 

The derived orbital decay for X— 7 is well within 
the range of values determined for other HMXBs 
like CenX-3 ((-1.738 ± 0.004) x 10~ 6 yr" 1 , Na- 
gase et al. 1992), SMC X-l ((-3.36±0.02) x 10~ 6 
yr"\ Lcvine et al. 1993) or LMC X-4 ((-9.8 ± 
0.7) x 10~ 7 yr"\ see e.g. Safi Harb et al. 1996; 
Levine et al. 2000). Such rapidly decreasing pe- 
riods in HMXBs are most likely caused by tidal 
interaction between the compact object and its 
massive companion. As the orbit decays the Roche 
lobe will descend into the companion's atmosphere 
and mass transfer will increase to super-Eddington 
rates over a relatively short time scale. In the end, 
the compact object is expected to spiral into the 
envelope of the companion and in this way termi- 
nate the high-mass XRB phase of the evolution 
(see e.g. Levine et al. 1993). 

4.2. The optical companion 

The HST WFPC2 images clearly resolve the 
dense OB association HS 13 (Humphreys & 
Sandage 1980). The optical counterpart is lo- 
cated to the North and is one of a pair of stars 
with similar luminosities: these stars are sepa- 
rated by ~0."9. While the suggested counterpart 
is presumed to be a single source, the source to 
the South is a blend of at least 2 stars (elon- 
gation from SSE to NNW). The HST observa- 
tions were carried out within one hour on October 
25, 1995 corresponding to binary phase 0.76 us- 
ing the ephemeris given in Sect. 4.1. Phase 0.76 
corresponds to the second maximum of the el- 
lipsoidal light curve. Assuming the parameters of 
the optical light curve by PMM2004 (sinusoidal fit 
with 0.033 mag amplitude) the correction to X-ray 
eclipse is +0.066 mag resulting in corrected mag- 
nitudes me of 18.9 mag, and my of 19.0 mag for 
the optical counterpart to X— 7. Magnitudes and 
colours are very close (within 0.1 mag) to values 



deduced by PMM2004. X-7 is also included in 
the UBVRI photometry of stars in M 31 and M 33 
from the survey of Local Group galaxies, recently 
published by Massey et al. (2006). They give a 
position (registered with the USNO B1.0 frame) 
of RA J200 o = 0l' i 33 m 34!18,(5j 20 oo = +30°32'll'/5, 
about 0'.'3 ESE from the position in our regis- 
tered HST image position (see Fig. 5). Their X— 7 
magnitude and colours (my = 19.087 ± 0.007mag, 
^(B-v) = -0.084 ±0.10mag, m (u _ B) = -1.057 ± 
O.llmag) coincide with those determined by us 
within 0.2 mag. Part of the discrepancy may be 
due to the sampling at different binary phases. 

The type and luminosity class of the star can 
be deduced from the absolute optical magnitude 
and colour during eclipse when we see the opti- 
cal surface that is mostly undisturbed by grav- 
itational effects, an expected accretion disk and 
heating by the X-ray source. To derive the ab- 
solute magnitude the measured brightness has to 
be corrected for the distance (-25.50 mag for the 
assumed distance of 795 kpc) and for interstel- 
lar extinction, the colour has to be corrected for 
reddening. These corrections have been estimated 
from the absorption of the X-ray spectrum in Sect. 

2.2. Based on the HST values the companion star 
should have an absolute My of —6.1 mag and 
(B — V)o of —0.3 mag. This corresponds to a giant 
star of spectral type 06III which has a tempera- 
ture of 39 500 K, a radius of 17 Rq and a mass 
well above 20 M© (see Appendix E in Carroll & 
Ostlie 1996). 

4.3. M 33 X-7, an eclipsing black hole 
HMXB 

With the new eclipse duration and the better 
determination of the companion type on the basis 
of the extinction combined with the colour excess 
corrections derived from the absorbing column of 
the X-ray spectrum of X— 7, we can significantly 
improve the mass estimate of the compact object 
compared to PMM2004. To do so we correct the 
B and V light curves given in PMM2004 for ex- 
tinction and colour excess. We then model these 
light curves using the "PHysics Of Eclipsing Bi- 
nariEs" program PHOEBE (Prsa & Zwitter 2005) 
built on top of the widely used WD program (Wil- 
son & Devinney 1971; Wilson 1979, 1990). By 
adjusting the semi-major axis of the binary sys- 
tem we kept the radius of the secondary star at 
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17i?0. We also fixed the temperature of the sec- 
ondary to 39 500 K. We then fitted the luminosity 
of the secondary and the mass ratio of the com- 
panions using the light curves in the B and V band 
simultaneously, assuming different inclination an- 
gles. We derive acceptable fits for inclination an- 
gles from 75° to 90°. Angles of 74° and smaller 
can be excluded as the companion star would over- 
fill its Roche lobe. The derived parameters are 
given in Table 3. The available optical data are 
not sufficient to discriminate between the allowed 
inclination angles. However, we expect that an 
inclination above 80° is more likely than an incli- 
nation in the range 75° - 80° where the compan- 
ion nearly fills its Roche lobe and unstable mass 
transfer would be expected based on the X-ray ob- 
servations. In the case of unstable mass transfer, 
the X-ray emission from X— 7 would most likely 
be much more variable on longer time scales than 
observed in all observations since those with the 
Einstein observatory. This implies a mass of the 
compact object in the system of greater than 9M Q 
and clearly indicates a black hole as the compact 
object in the system. 

As already discussed in PMM2004, further ar- 
guments for the black hole nature of the compact 
object in X— 7 come from the lack of X-ray pul- 
sations, the short term variability and the X-ray 
spectra. We discuss each of these properties in 
turn: 

(I) Pulsations are clear indicators of a neutron 
star as the compact object in a HMXB. In the 
power density spectrum analysis we did not detect 
significant periodic signals allowing a black hole as 
the compact object. However, this does not rule 
out the possibility that the compact object is a 
neutron star. Our unsuccessful Chandra periodic- 
ity search was limited towards short periods by the 
ACIS-I sampling time of 3.2 s. Also, with the sig- 
nificantly shorter sampling time of XMM.-Newton 
EPIC pn of 0.073 s, PMM2004 did not find signif- 
icant pulsations. 

(II) Low accretion rate sources should show 
a power density spectrum (PDS) with a broken 
power law (see Barnard ct al. 2005, and refer- 
ences therein). Above a luminosity of typically 
0.1 of the Eddington luminosity, the PDS should 
be flat. The short term fluctuations seen in the 
power spectral analysis of X— 7 are very small as 
one would expect for such a high accretion rate 



source. The unabsorbed luminosity of X— 7 in the 
0.3-10 keV band of >l.lxl0 38 erg s _1 in maxi- 
mum is consistent with a stellar mass black hole, 
but it does not constrain the mass. 

(Ill) HMXBs with a neutron star as the com- 
pact object normally show power law spectra 
(photon index 0.8-1.5) with a high-energy cut- 
off around 10-20 keV (see e.g. White et al. 1983; 
Makishima et al. 2000). Disk-blackbody spec- 
tra, on the other hand, suggest emission that is 
dominated by the inner accretion disk of a low 
mass XRB system, or - in case of a HMXB - 
the presence of an accretion disk surrounding a 
black hole emitting in the high state (e.g. Mak- 
ishima et al. 1986). As mentioned before, the 
variability of X— 7 outside eclipse most likely is 
caused by partial covering of the X-ray emission 
region due to material in the accretion stream or 
in the outer accretion disk. If we assume that dur- 
ing the brightest phase, the X— 7 disk-blackbody 
normalization corresponds to the innermost sta- 
ble circular orbit with radius r around a black 
hole, i.e. r = 3i? s = 9 x 1O 5 (M X /M ) cm (with 
R s = 2GM x /c 2 the Schwarzschild radius), we ob- 
tain M x > 2.4M . 

All of these results of the new Chandra obser- 
vations suggest that M 33 X— 7 is the first known 
eclipsing HMXB with a black hole as the compact 
object. 

Pooley & Rappaport (2005) advanced the idea 
of searching for eclipses in ultra-luminous X-ray 
sources (ULXs) to determine black hole masses 
and - most importantly - to separate intermediate 
mass from stellar mass black holes as the compact 
object in this exceptional class of X-ray sources. 
They proposed to compare the number of eclipses 
in the different kind of systems and predicted that 
more eclipses by far should be detected in stellar 
mass systems than in intermediate mass black hole 
systems. The orbital periods and other system pa- 
rameters would provide considerable insight as to 
the nature of the binary. While M33 X-7 is about 
a factor of five less luminous than the commonly 
accepted lower limit of ULX luminosities, it is still 
the first eclipsing black hole XRB system and - if 
in other observations detected in a factor of five 
brighter phase - may be the first stellar mass ULX 
candidate suitable for the proposed type of inves- 
tigations. 

As we discussed above, the resolved eclipse 
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ingress and egress in the Chandra energy band is 
not caused by the extent and the structure of the 
X-ray emitting region but by absorption in the 
density structure of the atmosphere of companion 
and the accretion stream. Watarai et al. (2005) 
investigated the asymmetry of black-hole accre- 
tion flows and suggested the presence of effects 
with time scales of seconds due to relativistic ef- 
fects close to the black hole. These should lead to 
asymmetries in the eclipsing light curves on time 
scales of seconds. The Chandra observations of the 
eclipse in- and egress of X— 7 show that such ef- 
fects are masked by much longer time scale effects 
discussed above. They may only be observable 
by sensitive observations at higher X-ray energies 
such as those where absorption effects in the com- 
panion atmosphere and the accretion flow vanish. 

5. Conclusions 

The ChASeM33 observations clearly resolved 
eclipse ingress and egress of the persistent eclips- 
ing HMXB M 33 X-7. Combining the eclipse 
egress times with Einstein, ROSAT, and XMM- 
Newton observations allowed us to improve on the 
orbital period which indicates orbital decay. The 
X-ray spectrum can best be described by a disk- 
blackbody spectrum with parameters that do not 
change significantly in different observations while 
the source luminosity changes by up to a factor of 
two. The short term variability of X— 7 can be 
described by a flat power density spectrum. No 
significant regular pulsations were found in the fre- 
quency range of 10" 4 -0.15 Hz. HST WFPC2 im- 
ages resolve the optical counterpart, which with 
the help of extinction and colour corrections de- 
rived from the X-ray absorption can be identified 
as an 06III star. The X-ray period and eclipse 
duration - together with fits to the light curve 
by PMM2004 of the optical companion imply a 
black hole with a mass above 9M Q as the most 
likely compact object in the system. M 33 X— 7 is 
the first eclipsing black hole HMXB. 

This research has made use of data obtained 
through the High Energy Astrophysics Science 
Archive Research Center Online Service, provided 
by the NASA/Goddard Space Flight Center, of 
the USNOFS Image and Catalogue Archive op- 
erated by the United States Naval Observatory, 



Flagstaff Station (http:/ /www. nofs.navy.mil/data/fchpix/), 
of data products from the Two Micron All Sky 
Survey, which is a joint project of the Univer- 
sity of Massachusetts and the Infrared Processing 
and Analysis Center/California Institute of Tech- 
nology, funded by the National Aeronautics and 
Space Administration and the National Science 
Foundation. The authors thank Michael Bauer 
for his introduction to the "PHysics Of Eclips- 
ing BinariEs" program PHOEBE. We thank the 
anonymous referee for his/her comments and sug- 
gestions for improving the manuscript. 
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Tabic 1: Chandra observations of the ChASeM33 



program covering M 33 X— 7. 



ChASeM33 


ObsID 


Obs.dates a 


Elapsed time 


ACIS 


Offax 


X— 7 binary 


Field 






(ks) 


CCD-ID 


(') 


phase 3. 


cycle b 


2 


6378 


2005-09-21 


112 


7 


17.9 


0.675-1.044 


-2,-1 


5 


7170 


2005-09-26 


41 


1 


1.6 


0.200-0.334 





5 


7171 


2005-09-29 


38 


1 


1.6 


0.041-0.165 


1 


5 


6384 


2005-10-01 


22 


1 


1.6 


0.832-0.905 


1 


6 


6386 


2005-10-31 


15 





7.2 


0.365-0.412 


10 


6 


7196 


2005-11-02 


23 





7.2 


0.928-1.002 


10,11 


6 


7197 


2005-11-03 


13 





7.2 


0.144-0.185 


11 


6 


7198 


2005-11-05 


22 





7.2 


0.798-0.868 


11 


6 


7199 


2005-11-06 


15 





7.2 


0.017-0.063 


12 


6 


7208 


2005-11-21 


12 


1 


7.3 


0.622-0.660 


16 


4 


6382 


2005-11-23 


73 


2 


9.7 


0.071-0.312 


17 


4 


7226 


2005-11-26 


25 


2 


9.7 


0.841-0.925 


17 



"date of start of observation 

6 with respect to eclipse center HJD 245 3639.119 and orbital period 3.453014 d (see text) 



Table 2: HST Observation Log of the M 33 X-7 
field. 



Dataset 


Date 


Time 


Exposure 

(s) 


Filter 


u2tr0301t 


1995-10-25 


13:40:16 


160 


F555W 


u2tr0302t 


1995-10-25 


13:45:17 


400 


F336W 


u2tr0303t 


1995-10-25 


13:54:17 


400 


F336W 


u2tr0304t 


1995-10-25 


14:04:17 


500 


F439W 


u2tr0305t 


1995-10-25 


14:15:17 


300 


F439W 



Table 3: M 33 X— 7 binary parameters for different 



inclinations determined with PHOEBE. 



Inclination 


SMA a 




M opt 


M x 


(degree) 


(Re) 


(M e ) 


(M e ) 




90 


38.1 


3.68 


49.1 


13.3 


85 


37.4 


3.88 


47.0 


12.1 


80 


35.8 


4.53 


42.4 


9.4 


75 


33.5 


5.78 


36.2 


6.2 



a semi-major axis of binary system 
6 mass ratio (Secondary over Primary) 
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M33 Field 4 

Bin time: 3.841 




10 3 0.01 0.1 

Frequency (Hz) 

Start Time 13697 13:86:34:393 Stop Time 13698 9:29:33:366 



Fig. 2. — Chandra ACIS power density spectrum 
of M 33 X-7 from ObsID 6382 using 23 intervals 
of 1024 3.241 s bins. 
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Fig. 4. — Deviation of eclipse egress times relative 
to best orbital period of 3.453014 d and Chandra 
eclipse egress ephemeris. The dotted lines give 
the shortest and longest periods allowed by the 
Chandra and ROSAT data. The parabola assumes 
a period of 3.45294 d at the Chandra epoch and 
orbital period decay rate of P or b/-Porb = — 4 x 10~ 6 
yr _1 (see text). 



M33 X-7 (outside eclipse) 




Channel energy (keV) 

Fig. 3.— Chandra ACIS spectra of M 33 X-7 
during high state. Data and the corresponding 
disk black body model are shown (see text). For a 
better representation the spectra were multiplied 
successively by factors of 1, 5, 5 2 , 5 3 and 5 4 from 
the bottom to the top. 




1:33:34.6 1:33:34.4 1:33:34.2 1:33:34.0 
Right Ascension (J2000) 

Fig. 5.— HST WFPC2 F439W 10"xl0" image of 
the HS 13 field with X-ray positions overlaid. The 
black circle gives Chandra position and error from 
Field 5. The red circle indicates the position given 
by PMH2004 based on just one radio SNR. Image 
registration uncertainty is indicated by the cross 
in the lower left corner. 
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M 33 X-7 X-ray light curve 



-I, 

u 



pi 

o 

o 



CO 



O 

o 



7208 



02 

CD 



CD 




7196 



7171 

7199 6382 
7197 



7170 




6386 



O - 



o 

CO 



CD 

o 



CD 



Hi 




ROSAT HRI 



0J 



1.2 



1.4 



Phase 



Fig. 1. — Light curve of the X-ray binary M 33 X— 7. (Above:) Chandra ACIS light curve in the 0.5-5.0 keV 
band. Individual observations are marked by their ObsID. (Below:) ROSAT HRI light curve (0.1-2.4 keV). 
Observations around eclipse on JD 2449571 are marked as red squares. 
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